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Abstract 



We derive a nonsymmetrized 8-band effective-mass Hamiltonian for quantum- 
dot heterostructures (QDHs) in Burt's envelope-function representation. The 
8x8 radial Hamiltonian and the boundary conditions for the Schrodinger 
equation are obtained for spherical QDHs. Boundary conditions for sym- 
metrized and nonsymmetrized radial Hamiltonians are compared with each 
other and with connection rules that are commonly used to match the wave 
functions found from the bulk kp Hamiltonians of two adjacent materi- 
als. Electron and hole energy spectra in three spherical QDHs: HgS/CdS, 
InAs/GaAs, and GaAs/AlAs are calculated as a function of the quantum dot 
radius within the approximate symmetrized and exact nonsymmetrized 8x8 
models. The parameters of dissymmetry are shown to influence the energy 
levels and the wave functions of an electron and a hole and, consequently, the 
energies of both intraband and interband transitions. 
PACS numbers: 73.20.Dx, 73.40.Kp, 73.40.Lq 

I. INTRODUCTION 

The 4x4 kp hole Hamiltonian for the wave function envelopes (so called effective- 
mass Hamiltonian), that takes into account mixing of the light- and heavy- hole bands, was 
obtained in Ref. p] using the perturbation theory. This multiband Hamiltonian has been 
employed for description of the hole states in bulk crystals! as well as in low-dimensional 
structures, in particular, in free-standing homogeneous quantum dots (QDs)IR The inclu- 
sion of the mixing with the snin-orbit split-off hole band leads to the 6x6 k p Hamiltonian 
which has also been appliedt™ to QDs. To consider the nonparabolicity of the electron dis- 
persion in narrow- and medium-gap semiconductors, it is necessary to take into account the 
coupling of the conduction and valence bands. Using the k p perturbation theory for bulk 
semiconductors with cubic lattice symmetry, the 8x8 k p model was developed in Ref. [7]. 
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This model explicitly includes eight bands around the T point of the Brillouin zone, namely, 
electron, heavy-, light-, and spin-orbit split-off hole bands (each of them is twice-degenerate 
due to the spin), and treats all other bands as remote. Along with more simple models, the 
8x8 kp Hamiltonian has been used to investigate different QDs (see, e. g. Refs. |8"Hl"Tl). 

Recently, one has begun to apply multiband effective-mass Hamiltonians to investigate 
elastic, electronic, and optical properties of multilayer nanostructures such as auantum- 
dot heterostructures (QDHs): CdS/Hg S0, InAs/GaAsSS, GaAs/AUGai^AsSEl, and 
CdS/HgS/CdS/H 2 0ll3tJ. However, it should be emphasized, that multiband kp Hamil- 
tonians are derived for homogeneous bulk materials, i.e. under the assumption that all 
effective- mass parameters are constant. This is important, because at a certain step of the 
derivation, wavenumbers k are declared as operators p/h that do not commute with the 
functions of coordinates. But, at the heterointerfaces of the multilayer nanostructures, there 
occurs an abrupt change of effective-mass parameters from their values in one material to 
those in the adjacent material. Inside a thin transitional layer that contains the heteroin- 
terface, the ordering of the differential operators and coordinate-dependent effective-mass 
parameters in the multiband Hamiltonian becomes crucial. In QDs with an infinitely high 
confining potential for electrons and holes, all components of the wave function vanish at 
the heterointerface, and there remains a possibility of applying the bulk multiband k p 
Hamiltonian straig htforwardlyihlHill There are two ways to proceed from QDs to QDHs. 

(i) The first way is to use an appropriate bulk multiband Hamiltonian for each con- 
stituent material separately, and then to match the obtained homogeneous solutions at the 
abrupt hetero junctions applying the connection rules (CRs) that are usually obtained by 
imposing the continuity of the wave function envelopes and of the normal to the heteroint- 
erface component of the velocity.§@ It should be underlined that this way is heuristic and 
nonunique. In Ref. [TBI the general CRs, that even do not require the continuity of the wave 
function envelopes, have been proposed for planar heterostructures. 

(ii) The second way (cf. Refs. |TD|-[n|) is to derive a multiband Hamiltonian valid for 
the entire heterostructure, including the heterointerfaces, and then, if material parameters 
change abruptly at some interfaces, to find the boundary conditions (BCs) for the solu- 
tions of the envelope function equation. To find these BCs, one should use the multiband 
envelope function equation (h — E^j ^ = at any point of the heterostructure, including 
the heterointerfaces, and integrate this equation over the volume of an infinitely thin layer, 
which includes the considered heterointerface. Thus, the BCs are derived starting from the 
requirement of continuity of the components of the wave function at the heterointerface. 

One can always choose the CRs physically equivalent to the BCs.il The both approaches 
(i) and (ii) are usually used when the wave function inside each layer of a heterostructure 
can be found analytically, for example in planar or spherical heterostructures. In case of 
an arbitrary shape of the heterointerface, the approach (ii) can still be used because, when 
the Hamiltonian is known for the entire heterostructure, one can find an overall numerical 
solution of the Schrodinger equation. 

A commonly used heuristic method to obtain a multiband effective-mass Hamilto- 
nian for heterostructures uses symmetrizationliHEl of the corresponding k p Hamilto- 
nian. This method consists in the symmetrical arrangement of the components of the 
momentum operator, that ensures the hermicity of the resulting Hamiltonian. Namely 
(3 p — > (f3(r) p + p f3(r)) /2 and (3f>iPj — > (pi /3(r)pj + pj (3(r) pi) /2, where /3(r) is a spatially 
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varying effective-mass or other material-dependent parameter which is usually considered a 
piece- wise constant, because in each layer of a heterostructure it has the value for a corre- 
sponding bulk material. The symmetrization has been applied to QDHs in Refs. [I2| , [i3| . |T7 . 
An essential fault of the symmetrization is that it is not a necessary condition for the 
multiband Hamiltonian to be hermitian. Besides that, as will be seen below, some intrinsic 
properties of the heterointerface, such as reducing the symmetry of the problem and smooth- 
ing the abrupt change of the effective-mass parameters at a heteroj unction, are completely 
neglected in the symmetrized Hamiltonian. 

Burt has derived^UHll the exact envelope-function equations for a heterostructure. The 
order of the components of the momentum operator arises as a part of that derivation. This 
theory has been used by Foreman to explicitly write the 6x6 (Ref. |28|) and 8x8 (Ref. |29l) 
effective-mass Hamiltonians for planar heterostructures. General rules for constructing the 
valence-band effective-mass Hamiltonians with a correct operator ordering have been de- 
scribed in Ref. [JU] for the heterostructures with arbitrary crystallographic orientations. In 
Ref. 31,32], correct boundary conditions for planar heterostructures with wurtzite symmetry 
have been presented. Comparing the conduction- and valence-subband dispersion of a planar 
quantum well, calculated using the BCs following from the exact nonsymmetrized and from 
the symmetrized effective-mass Hamiltonians, it has been shown that the former BCs 
physically reasonable results, while the latter BCs can produce nonphysical solutions.1 
More recently, these two sets of BCs for a planar quantum well have been examined within 
the tight-binding approach.EI The result of the comparison allowed to give preference to the 
nonsymmetrized model resulting from Burt's derivation of the envelope-function Hamilto- 
nian, which was shown to give reliable results even when the well and barrier effective-mass 
parameters were very dissimilar. 

In the present paper the envelope-function representation of Refs. P^-f2T| is used to con- 
struct the nonsymmetrized 8-band Hamiltonian for an arbitrary 3-dimensional heterostruc- 
ture. As an application, the electronic structure of two-layer HgS/CdS, InAs/GaAs, and 
GaAs/AlAs spherical QDHs is investigated as a function of the dot radius. 

It should be mentioned that the spurious solutions ( "oscillating" @ states and "gap" 
states0H) did not become apparent in the aforementioned QDHs. However, such solutions 
may appear for a different set of parameters. 

The results of the calculation are compared with those obtained from the symmetrized 
8x8 Hamiltonian. The rest of the paper is organized as follows. In Sec. II the derivation 
of the nonsymmetrized 8-band Hamiltonian for a QDH is presented. The corresponding 
radial Hamiltonian for a spherical QDH is obtained in Sec. III. In Sec. IV the BCs for both 
symmetrized and nonsymmetrized radial Hamiltonians are compared with each other and 
with commonly used CRs. The results of the numerical calculation for spherical QDHs 
are obtained and discussed in Sec. V. Conclusions are given in Sec. VI. The 2x2 electron 
and 6x6 hole energy- dependent nonsymmetrized Hamiltonians for a QDH, as well as radial 
Hamiltonians and corresponding BCs for a spherical QDH are found in Appendix B from 
the nonsymmetrized 8-band Hamiltonians. 
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II. NONSYMMETRIZED 8-BAND HAMILTONIAN 



We begin our derivation with the nonsymmetrized 8-band effective-mass Hamiltonian for 
a heterostructure, when the spin-orbit coupling is "turned-off " . In the Bloch function basis 
\S), \X), \Y), \Z) this Hamiltonian is represented in the following forml! 



Ha 



2m 



e c + kak - [yik x + k 



rV 2 



V-2 



% - (yxk z + k, 



V2 



v 2 k x + k x v^j e' v - k x Pik x - k^p h k^ -3 (k x lt K + h^kx) -3 (^7^ h + k z ^k x ^ 
o y v 2 k y + kyV^j -3 (k x jsk y + k y j£ k x ^ e' v - kyfyky - kyf3 h ky -3 (^7^ & z + hl^h) 
(v 2 kz + k z vi) -3 {k x ^k z + k z j£ k x ^j -3 {k y ^k z + k z j£ k y ^j e' v - k z (3ik z - k z /3 h k z 

(1) 

where e' v = e v - 5/3, k = -iV, k^z = k- k x , y , z , A = 7i + 472, Ph = li ~ 2l 2 , 

vi = v + t, v 2 = v-£, (2) 



73 = 7s + X, 7 3 = 7s - X- 



(3) 



£ and x are called the dissymmetry parameters, because when £ = and x = 0, the 
Hamiltonian ([I]) becomes symmetrical. The explicit form of the parameter £ = {y\ — v 2 )/2 
follows from the formulae (|Sl|) and ([A 2D in Appendix A, and the parameter x is determined 
byEq. 

When the spin-orbit coupling is "turned-on" , the considered 8-band Hamiltonian is rep- 
resented in the Bloch function basis \S]}, |Xt), |FT), \Z]}, \Si), \Xi), \Y[), \Z[) as 



H A 
H 4 



H x 



(4) 



where Ha is defined by Eq. ([]]) and the spin-orbit Hamiltonian H so has the forme 
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(5) 



In Eq. (P, mo is the free-electron mass, -E c = h 2 e c /2mo is the energy of the conduction 
band (CB) minimum, E v = h 2 e v /2mo is the energy of the valence band (VB) maximum, 
A = h 2 5/2mo is the spin-orbit splitting of the VB, and V = hv/2mo is the Kane velocity 
(V = —iK(S\k z \Z)/mo). Contributions of remote bands to the hole effective masses are 
written in terms of the "modified" Luttinger parameters 7i = 7i — E p /3E g , 72 = 7^ — 
E p /6E g , and 73 = 7^ — E p /6E g , where E g = E c — E v is the energy gap, E p = 2m V 2 is the 
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Kane energy, and 7, L (i = 1,2,3) are the Luttinger parameters of the VB. Parameter a can 
be evaluated through the experimentally determined CB-mass m c using the relation 
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2 1 

E~ 9 + E g + A, 



(6) 



It is worth noting that all parameters entering the Hamiltonian (Q) are coordinate-dependent. 
In a heterostructure these parameters abruptly change from their values in one material to 
the corresponding values in the adjacent material, therefore they are piecewise-constant 
functions of r. Although not symmetrical, the Hamiltonian H s is hermitian as seen from 
Eq. (|J). The parameters £ and \ (see Eqs. (0) and (0)) are responsible for the nonsymmet- 
rical form of the Hamiltonian (|j). The symmetrized 8-band Hamiltonian can be obtained, 
therefore, taking £ = and x = 0- 

In Eq. (|]), 73" is the contribution to 73 from the r x and T 12 remote bands, while 73" is 
the contribution to 73 from the T 15 and T 2 5 remote bandsS Neglecting small contributions 
from the r 2 5 remote bands, the parameter x(r) is determined as 
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X (r) = (2 72 (r) + 3 73 (r) - 7l (r) - l)/3, (7) 

i.e. it is explicitly defined by the effective-mass parameters of the bulk model. It is seen 
from Eq. (|7j) that in a homogeneous medium when 7 j(r) are constants, x( r ) is also a con- 
stant, and therefore, it cancels from Eq. (p. Consequently, x( r ) is a specific function of a 
heterostructure, which gives a nonzero contribution to the Hamiltonian only at the heteroin- 
terfaces. The value of this contribution at the point ro of a heterointerface is proportional 
to x( r o + e n) — x( r o — e n), where e n is an infinitesimally small vector, normal to the het- 
erointerface at the point r . Parameters 17 (r) and f 2 (r) of the Hamiltonian_i|l|) , which can 
be obtained from the general effective mass equations as derived by Burtt3, are given in 
Appendix A. In the definition (|2]), the functions ui(r) and ^2(1") are subdivided into the 
symmetric v(r) and antisymmetric £(r) parts, where £(r), like x( r ) above, is a specific pa- 
rameter of a heterostructure. In general, £(r) is a piecewise-constant function of r. The 
necessary and sufficient condition for £(r) to give a nonzero contribution to the 8x8 Hamil- 
tonian only at the heterointerfaces, and to become a constant in the homogeneous medium, 
simultaneously with v(r), is 

£(t) = ccv(t), (8) 

where the coefficient of proportionality is constant over the entire heterostructure. Eq. (§) 
is the general form of £(r) only for a two-layer heterostructure. For an iV-layer heterostruc- 
ture, there can be N — 1 independent constants - one for each heterointerface. Each constant 
for a given heterointerface can be found experimentally considering a two-layer heterostruc- 
ture (see Sec. V). 

In order to diagonalize the spin-orbit Hamiltonian H so , it is convenient to carry out a 
unitary transformation of the Bloch function basis 15*1), |A|), |F|)> \Z^), |A|), \Yl), 
\Z[) into the following Bloch function basis! 

n< L/2,l/2 = l*ST/> 

(9) 

M l/2-l/2 — l&l)> 
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"3/2,3/2 ~~ 


-L(|xt) + i|YT>), 


"3/2,1/2 — 


±Mx\) + i\Y\)-2\Z\)) 
V6 v 1 


V 

3/2,-1/2 ~~ 


^(ixt>-<m>+2|zi>) 


V 

3/2,-3/2 


' (ixi)-zin)), 


"1/2,1/2 ~~ 


-L(l^>+iin> + |2T>), 


V 

1/2,-1/2 - 


^(ixt)-.in)-i^)), 



(10) 

Un 



u 

where itj and are the Bloch functions of the conduction and valence bands, J is the 
Bloch function angular momentum, and /i = J z is its ^-component. The 8-band Hamiltonian 
H in the new basis can be obtained now after a unitary transformation 

H = U*H 8 U T , (11) 

where U is the matrix of the transformation from the basis {|5"T), |^t), |Y|), \Z^), \Si), 

1^1), l^l}} tO the new basis {Uy 2 ,l/2i M l/2,-l/2' M 3/2,3/2' M 3/2,l/2> U 3/2,-l/2' M 3/2,-3/2' 

u i/2i/2' ""1/2-1/2}- ^ * s defined by Eqs. (Q) and ([X0|) . Performing the transformation (^Tj), 
one obtains: 

£ C + T iV, ,/|v„ -=l/_! -=v„ ,/=v_, 

+ T if-V -v., i^/fvi ±Vi 

-iV? e v -P-Q -S -R ^=5 zV2P 

v2 

^2 y|ltf ^ -5t e.-P + g -C -P iv^Q -z^£ 

2mn 



-yi : -i?t S* e v -P*-Q* iy/2B) -^S* 

v2 

^Lvj -i^fv? ^ -iV2Q -z t /fs r -ij2R <-P C 

^1 -iV2tf zyfst -zV2Q* ^ C* e»-P* 

(12) 



7^ 





#° 


-V-x 


-P 





-c 


-P 


-P* + Q* 


S T 


S* £ v 


- P* - Q* 




-iV2R 


-zV2~g* 


\/2 



where e" = e v — 5, 



v i = 2 ( v ^+ + ^2) , Vli = - (yxk- + A;_f 2 ) , 
^0 = g + ^2) , T = k + ak_ + k^ak + + k z ak z , 

P = fc + ( 7l - 2 X )fc- + fc-(7i + 2x)^+ + ^7i^, 
Q = k+{j2 ~ x)k- + M72 + - 2^72^, 
i? = v/3 (A; + (72 - 7 3 )A; + + k-(j 2 + 73)^-) , 

S 1 = -iy/Q (L(7 3 + x)A) 2 + M73 - > 

e = -iVe (L( 7 3 - |A + ^(73 + |)L) , 

C = -z 2 V2 (k-xk ~ kxk-) ■ (13) 

In Eq. (|12"D , daggers (f) denote the hermitian conjugation, i.e. = (^4 T ) (it is impor- 
tant to note, that (vi^k^j = k T v\^ and (vi^k^ = k z V\^). Unlike the bulk 8x8 kp 
Hamiltonian^ where the matrix element C is zero, in the Hamiltonian fll2"|) the Bloch func- 
tions W3 , 2 y 2 an d ^3/2 -1/2 are coupled with each other, as well as the functions u\, 2 u 2 and 
u\, 2 _ 1 / 2 . As seen from Eq. ([13]) , this coupling arises because of the dissymmetry parameter 
X, which can reduce, in this way, the symmetry of the problem. The Hamiltonian so ob- 
tained can be used to investigate electronic properties of quantum-well, quantum-wire, and 
quantum-dot heterostructures. 



III. 8-BAND HAMILTONIAN FOR A SPHERICAL QDH 

To study the electronic structure of spherical QDHs, the spherical approximation! (i.e. 
I2 = 73 = 7 L ) can b e applied. If we take 7 L = (27^ + 373 )/5, then the quantum states 
so obtained are correct to the first order of the perturbation theory. Using the relations 
between the Luttinger parameters 7^ and the "modified" Luttinger parameters 7$ we have 

7 = (272 + 3 73 )/5, (14) 

and according to Eq. ([?]) 

X=(5 7 -7i-l)/3. (15) 

In spherical QDHs, where all effective-mass parameters depend only on the absolute value 
r of the radius-vector, electron and hole states are eigenf unctions of the total angular mo- 
mentum j and its ^-component m = j z . Therefore, the electron or hole wave function can 
be written as a linear expansion in the eight Bloch functions Uj . 

*i,m(r) = E «^ + £ ^T(r) «V ( 16 ) 
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where the envelope functions Fj^/'^' m (r) are defined in the chosen Bloch function basis 
(0) as 



F 1/2™( T ) = E ^l/2, K l,X R y2,l( r ) ^,a(0> <t>)i 
l,X 

=« M/2 E^,A%MnA(M), (17) 

l,X 

= ^ E C^ w R^(r) Y ljX (9, 0). 
i,x 

Here, R c jf' J (r) are the radial envelope functions, Cj'"! ; A are the Clebsch-Gordan coeffi- 
cients, and Yi t \(9,(j)) are the spherical harmonics. Noting that in the matrix represen- 
tation of the Hamiltonian (m) u\j 2 y 2 — (l , . . . , Uy 2 -1/2 = 
(0 l) one can rearrange Eq. ([16]) into the form 

i+1/2 

«W*) = E ^M^CM) 

2=5-1/2 
j+3/2 

+ E ^(r)3T^,^) (18) 
1=3-3/2 

3+1/2 

+ e KU r )y$sQ>®> 
1=3-1/2 

where the 8x8 matrices yji'^' m (9,(f>) next to the eight radial envelope functions (r), 



for a given j, can be found by comparing Eq. (TS) with Eqs. ( JIB] ) and (0). Now, integrating 
over the angular variables 9 and 0, it is possible to obtain the radial Hamiltonian 

% = I (yj>T& y h u"\0- <P) dn, (19) 



'J,l 

corresponding to the radial Schrodinger equation 

E^ R%{r) = E 3 R%(r\ (20) 

b,J,l 

where Ej is the electron or hole eigenenergy to be determined, b = c or v . The Hamiltonian 
([19|) does not depend on m, because within the spherical approximation the energy spectrum 
is degenerate with respect to the ^-component of the total momentum. 

After some algebra, we derive the following relations for the spherical harmonics: 

k + Y l>x {9, 0) = Cfti^.! B+ Y l+1)X+1 (9, <P) 

+ Cz-l,A+l;l -1 B l ^-l,A+l(0, 0), 

kJY hX {9,cj>) = -^; A liA _ 1;1>1J B+^ +liA _ 1 (0,0) 

-Cj-l,A-l;l,l ^ *i-l,A-l(0, 0), 
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KY, 



Cf+i,A ; i,o Yi+i, 



+ Q-1,A;1,0 ^ ^-1, 



(21) 



where 



1 + 1 



21 + 1 



l 5 



b; 



i 



21 + 1 



.4 



(-i) 



A 



(p) 



9 l + l/2-p/2 



or 



(22) 



Using the relations (pi ), the radial Hamiltonian fllTf ) can be obtained in an explicit form. 
If we choose the following order of the radial functions: Ru2j-i/2i -^3/2j+i/2' ^-3/2j-3/2> 
^L/2j+i/2> i? V2,j+i/2> R 3/2,j-i/2, i2 3/2j + 3/2> R T/2,j~i/2> then the 8x8 matrix of the Hamilto- 
nian 7ij takes the form 



7Yo 



(i) 



o n 



o 

(-i) 



(23) 



Here, TCj is the 4x4 Hamiltonian of the "even" states and TCj is the 4x4 Hamiltonian 
of the "odd" states. It is seen, that the parity p (p — 1 for even states and p — — 1 
for odd states) is conserved in the spherical approximation, even when the Hamiltonian is 



not symmetrized. The obtained radial Hamiltonian for the radial functions R\ 



/2J-P/2' 



K. 



v,3 



R 



2 ,/+/, 2' -"-3/2,j-3p/2> an d #l/2,j+p/ 2 naS tne ^° rm 
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(p) 



& 2 



2m r 



z i-p/2 
p »(p) 

a j °j-p/2 

b p 

\v^ P /2 



P y((-p) 



fc« i- / j+p / 2 



a 



C j ^-j+p/2 
(-P) 



^p jCp) 
p i^7->(p) 

S °j 'S'-3p/2 



J "J ' ~j+p/2 



_i_ -p(p) , jp qCp) 

u ' i-3p/2 ~ c j ^j-3p/2 

pV2 6f<2 3p/2 



Pvi *^i+p/2 
pV2b p j TZ H ' ] 



p n( - p) 
i+p/2 



i+p/2 



^- J + ^p/2/ 



(24) 



where = + 3r)j, = V 3 ( 1 ~ °i = 1 ~ 3r ?j > = ( 2 -? + 1 ~ ?)» 



(p) 



(25) 



ft 



(p) 



^z+p 7 A i ■ 



Introducing the operator 



Af\p) = -A\-;>pA 



(-P) a a(p) 



(26) 
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we can represent %, V l ■ , and Qj as 

^_ (/ + l)Ai 1) (^) + /Ai- 1) (g) 
' 2/ + 1 

^, = (^i^-^+z^-ax) + Afp)(2x); (27) 

W _ g ~ l/2)Aj 1) ( 7 ~ X) + (j + 3/2)Aj- 1) ( 7 - x ) (p) 

«i - 2TTI ' {x) ' 



Inside each spherical layer, the radial Hamiltonian for a spherical QDH coincides with 
the bulk radial Hamiltonian from Ref. § for a spherical QD, when the following denotations 



for the radial functions are used: 








E>c;i _ 
"l/2,7-l/2 " 


TD+ 


•"1/2^+1/2 


p— 

"-c,7 ' 


r>v;j _ 
-"-3/2J+1/2 


£>+ 


r>v;j 

"*3/2,i-l/2 


= *^M,7' 


p v J _ 
- rL 3/2,i-3/2 ~~ 


P+ 


pv;j 

"-3/2,7+3/2 


= 


"l/2j+l/2 - 


p+ 


pf;j 

"-1/2.J-1/2 


p— 



(28) 



Therefore, in order to find the radial wave functions inside each spherical layer, one 

can use the same technique as in Ref. ^| When the wave functions inside each spherical 
layer are known, the BCs should be applied to match the wave functions from two adjacent 
layers. 



IV. BOUNDARY CONDITIONS FOR A SPHERICAL QDH 



When considering the multiband models for planar heterostructures, the BCs for the wave 
function are often obtained by integrating the Schrodinger equation across the heterointer- 
face and assuming the continuity of the wave-function envelopes. SHIll. The resulting BCs 
are of the following form: 



B\ 



hO - 



B 



=+0 



(29) 



where A and B are two materials separated by the heterointerface z — 0, and j z is the normal 
to the heterointerface component of the current operator. The aforementioned integration is 
actually justified only for Burt's envelope-function equations, because only these have been 
shown to be valid at the heterointerface. Analogously to the case of planar heterostructures, 
for spherical QDHs one integrates the radial Schrodinger equation 



'nf - e) 



(p)\ p(p) 



lh 3 



0, R 



(p) 



/ n l/2j-p/2 \ 
p«!7 

- a 3/2j'+p/2 

pv;j 

"-3/2,i-3p/2 
\ U l/2,j+p/2 / 



(30) 
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where Tij is defined by Eq. (24) and Ej P> is the eigenenergy. This integration is carried 
out across the point r = a, where r = a is the spherical heterointerface that separates two 
materials: A (at r < a) and B (at r > a). Including the continuity of the radial wave 
function Itf , the required BCs have the form: 



n(p) 



.4 

' I A 



r=a—0 



r=a—0 



R 



(p) 



r=a+0 



B 



(31) 



r=a+0 



Here the radial component of the current operator Cf^ is obtained from the radial Hamil- 

tonian TCj using the following procedure. In those terms of the Hamiltonian (|24]) that 

contain the operator A^f\ the utmost left-hand-side A^ is replaced by — p (in conformity 
with Eq. (|22"D); the rest of the terms are set to zero; the result is multiplied by 2i/h. Thus 
we find 



ih 

m 



-a: 



d_ 

Or 



P 



A {v - o 



-p 



q (v - o 



-p 



2^6 3 2V6 



27f ( ^° 



'- % (v + -pa^ 3l A%) 3 7i I + cf 7 D r - 3ft f 



2^ 



py/2aShD r + fi 



ej,p 



X 



yftqiA) 



(p) 



-3p/2 



•%/2&?7^ P) 



J i j+p/2 

-i:7- + 2 /3- 



<9r 



r 
(32) 



where 



(33) 



While the radial Hamiltonian is hermitian, the radial component of the current operator 

is not (as seen from Eq. fl32|) ). 
It is important to compare the obtained BCs (|3~ID, ( |52"D with the commonly used CRs 
that the wave function and the normal component of the velocity are continuous at the 
heterointerface0Ei. The velocity operator V is defined as 



h 



H, r 



1 dH 



(34) 



where the Hamiltonian H has been determined earlier by Eq. (112]) . Therefore, the normal 
component of the velocity operator is obtained as follows: 



V r = ~- 



1 r dH 1 



h r dk 



- n 4 



dH 

dk 



+ n_ 



dH 

dk. 



+ n z 



oh" 

dk Zl 



(35) 



where 
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x + iy x — iy z 
n + = TpT* n - = TpT* n z = -. 36 

The differentiation of the Hamiltonian can be realized in such a way that 

Vr = Vr + K R , (37) 

where 

~ L <m 1 ( d H L{R) dH L ™ dH L W\ , , 

K ( } = r% + ™- + n z — • 38 

K V <9A; + <9A;_ dk z ) 

Here (-ff R ) denotes the Hamiltonian, in which the right (left) -hand operators k + , 
and fc 2 are treated as c- numbers, i.e. only the left(right)-hand operators k + , k_, and k z 
are differentiated. Using the explicit form of the Hamiltonian H (see Eq. (|12D ) , one finds 
that V r L can be obtained multiplying by 1/h the Hamiltonian H, in which all the left-hand 
operators k + , k_, and k z are replaced by n + , n— , and n z , correspondingly, while all the 
terms that do not contain the former operators are set to zero. It can be also shown that 

K R (^x) = rV r L (-^-x), (39) 

where the operator f draws all effective- mass parameters through the operators k + , and 
k z to the utmost right-hand positions. 

The radial velocity V^' L is obtained from V r L in the same way as the radial Hamiltonian 

Hj was obtained from H, i.e. by the definition fllTt ). For n + , n— , and n z (see Eq. fl36])) 

the expressions similar to ( |2T|) are valid if one replaces —i Af' by —p. Therefore, Vj L can 

be found by multiplying the Hamiltonian nf by i/h, replacing all the left-hand operators 

by — p, and setting all the terms, which do not contain the operator A^f\ to zero. This 
procedure results in 

Vf' L = (40) 

where Jj has been defined by Eq. (|32|) . Using Eqs. (|37|), pD]), and (|39"D one obtains 

)(p) _ I / -fteVt ,^ _l -r(p) 



Vf = -(J, W (e,x)+r^(-^,-x)j. (41) 

Considering the explicit form of the matrix J~j we see that all the terms, containing pa- 
rameters £ and Xi responsible for the nonsymmetrical form of the Hamiltonian, cancel. 
Consequently, we obtain 

y p) = 1 -^jj P \t = 0,x = 0). (42) 

If all the effective-mass parameters are piecewise-constant functions of r, then at any point 
of the heterostructure except the spherical heterointerfaces, we can replace f by 1, and 
therefore 
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v 3 



(p) 



j 7 (p) (e = o,x = o). 



(43) 

It is clearly seen now, that the commonly used CRs for spherical QDHs§@ are the same as 
BCs (|3~1~D, (|32|) obtained from the symmetrized Hamiltonian (£ = 0, % = 0). Like , the 
radial component of the velocity operator fl4"2|) is not hermitian when f = 1. Therefore, in 
order to prove that both current and velocity are conserved simultaneously, we should verify 
whether the real parts of the current density and of the velocity density are the same, i.e. 
we should check whether the equality 



Rc 



(P) R (P) 
J ^3 K 3 



Re 



(44) 



holds true. Here Rf is the radial wave function defined by Eq. (|30|). Substituting Eq. (|32|) 
into the left-hand-side part of Eq. (|44D , we see that all the terms containing the parame- 
ters £ and x cancel, because their contribution to the current density is purely imaginary. 
Therefore, in conformity with Eq. fl43|) , the equation (|44]) is proven to be fair. 



V. RESULTS OF THE CALCULATION AND DISCUSSION 

In this section we investigate the electronic structure of three spherical QDHs with dif- 
ferent values of the energy gaps: a zero-gap semiconductor embedded into a wide-gap semi- 
conductor (HgS/CdS), a narrow-gap semiconductor embedded into a medium-gap semicon- 
ductor (InAs/GaAs), and a medium-gap semiconductor embedded into a wide-gap semicon- 
ductor (GaAs/AlAs). Note that in these widely used experimentally relevant materials, the 
effective- mass parameters are substantially different. The bulk 8-band parameters of the 
used III-V and IV- VI materials are listed in Tables | and [TI|, correspondingly. For electron 
and hole levels, obtained within the spherical 8-band model, we use a common notation: 
nQ^ denotes an electron state and nQj denotes a hole state, where n is the number of the 
level with a given symmetry and Q = S, P, D, . . . denotes the lowest value of the momentum 
/ in the spherical harmonics of Eq. QT%D in front of the CB Bloch functions for an electron 
state and in front of the VB Bloch functions for a hole state, i.e. Q = j — p/2 for an electron 
and Q = min(j + p/2, \ j — 3p/2|) for a hole. 



A. Electron energy levels 

The electron energy levels of the HgS/CdS, InAs/GaAs, and GaAs/AlAs QDHs are 
depicted in Figs. correspondingly, as a function of the quantum dot radius a. The value 
of the spin-orbit splitting of electron energy levels is small (~ 3 meV) for all considered 
QDHs (see Table |T| and Appendix B). Therefore, only the lowest level of the pair of split 
levels is shown in Figs. [3]-|3]. For all examined QDHs, the lowest level of such a pair is the 
level with the least total momentum j. 

Analyzing Figs. [1|-|3] we arrive at the following empirical formula, which determines the 
energy shift of all electron levels with n — 1 when a nonzero value of the parameter \ is 
considered: 
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E e -E e \ x=0 = 2&1-X2). (45) 

x mo or 

Here the indices "1" and "2" denote interior and exterior materials, correspondingly. For 
n > 1 the value of this shift is much smaller, than for n — 1. It is seen that the shift 
E e — E e \ is about 16 meV for HgS/CdS (a = 2 nm) and about 4 meV for InAs/GaAs 
(a = 4 nm) and GaAs/AlAs (a = 3 nm) QDHs (see Table [TTI|). As provided by Eq. (f45|), 
the value of this shift is inversely proportional to the square of the quantum dot radius. 
Consequently, for large QDHs one can use with high accuracy the symmetrized with respect 
to x Hamiltonian, and the expression (|45| ) is the measure of accuracy. If Xi < X2, then the 
nonsymmetrized energy level lies higher than the symmetrized one (see Figs. [I], ^), and if 
Xi > X2, then the nonsymmetrized energy level lies lower than the symmetrized one (see 
Fig-ID- 

Gray bands in Figs. reflect the change of the parameter (see Eq. (|J)) from 1 to 
— 1. The chosen interval includes the following specific values of eg: eg = (symmetrized 
Hamiltonian) and eg = ±1 (see Appendixes A and B). With such a change of c^, the electron 
energy increases in Figs. [I], § and decreases in Fig. [|. Therefore, the shift of an energy level 
with n = 1 with respect to the level position when = can be estimated by the formula: 

E e -E e \ ee=0 = bc i ( y /E^- 1 /E^) (6 > 0), (46) 

where we take into account Eq. (|8|) and the fact that v\$ ~ J 1 E p ^^- E p is the Kane energy 
(see Tables | and (TI|). For n > 1 the shift E e — E e \ becomes much smaller, than for n = 1. 
The parameter b in Eq. (f46|) decreases with increasing a and with increasing the energy gap 
in the interior material. Such a behaviour of the parameter b is connected with the fact 
that it is proportional to the value of hole radial components of the electron wave function 
at the heterointerface. It is clear now that the observed strong dependence of the energy 
levels in the HgS/CdS QDH on (see Fig. [TJ) is due to the large value of the difference 
\jE Pi i — J 1 E p 2 and to the zero energy gap in HgS. For two other QDHs the dependence of 
the energy levels on is a few times weaker, than that for the HgS/CdS QDH. 

Hence, Eqs. (|^) and (|4"6| ) allow us to estimate corrections to the eigenenergies due to 
the replacement of the heuristic symmetrized Hamiltonian with the nonsymmetrized Hamil- 
tonian avoiding complicated calculations. It follows from Eq. (^) that such corrections rise 
with decreasing the quantum dot radius as 1/a 2 . Therefore, one should use the nonsym- 
metrized Hamiltonian for description of quantum dots with small radii. 

B. Hole energy levels 

All the hole energy levels of 5*- and P-types in the HgS/CdS QDH, with j = 3/2 in 
the InAs/GaAs QDHs and with j = 3/2 in the GaAs/AlAs QDHs are depicted in Figs. £|, 
|5] and ||, correspondingly, as a function of the quantum dot radius. It is seen from these 
figures that the empirical formula (^) holds for hole levels, too, i.e. for a nonzero value of 
X, the hole energy levels shift in the same direction as the electron energy levels do. For the 
HgS/CdS QDH (a = 2 nm) the shift of the hole ground state level is about 9 meV, what is 
smaller than the shift of the electron ground state level. At the same time, for InAs/GaAs 
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(a = 4 nm) and GaAs/AlAs (a = 3 nm) QDHs the shift for the hole ground state level is 
almost the same as that for the electron ground state level (see Table |T|). For the higher 
hole levels (n > 1) the value of the shift under consideration decreases with increasing n 
much weaker than it does for the electron levels. 

The dependence of the hole levels on the parameter eg is substantially different from such 
a dependence for the electron levels. The formula ([46]) can be approximately applied here 
only for the level 1S^\, which is the hole ground state energy for all examined QDHs. It 
is seen that this energy level strongly depends on even for InAs/GaAs and GaAs/AlAs 
QDHs. All the other hole energy levels under analysis depend on very weakly, and such 
a dependence is revealed only in Fig. [| for the HgS/CdS QDH. 

C. Electron and hole wave functions and pair energies 

In Figs. |7|, ||, and |9], the S'-components of the radial wave functions of the electron 
ground state (lSy 2 ) and of the hole ground state (lS^j\) are depicted for HgS/CdS (a = 

2 nm), InAs/GaAs (a = 4 nm), and GaAs/AlAs (a = 3 nm) QDHs, correspondingly. It 
is seen that in all these QDHs, the hole density in the interior material is higher than the 
electron density, and the electron density is higher in the exterior material. It is also seen 
that when changes from 1 to —1, the electron density in the centers of the HgS/CdS and 
InAs/GaAs QDHs increases and the hole density decreases. The opposite trends of behavior 
of the electron and hole densities are observed in the center of the GaAs/AlAs QDH. The 
abrupt change of the derivative of the electron radial component with the change of is 
well seen at the heterointerfaces of all QDHs under consideration. At the same time, the 
derivative of the hole radial component changes smoothly. The contribution of the hole 
radial components to the density of the electron state (at c^=0) is as high as 33 % for 
HgS/CdS, 20 % for InAs/GaAs, and 14 % for GaAs/AlAs QDH. Such contributions show 
that the nonparabolicity of the electron dispersion law is substantial even for the QDs of the 
medium-gap semiconductors (GaAs) and certainly should be taken into consideration when 
the QDs of the narrow-gap semiconductors (InAs) are investigated. The contribution of the 
electron radial component to the density of the hole state (at cg=0) is 6 % for HgS/CdS, 
1 % for InAs/GaAs, and 1 % for GaAs/ Al As QDH. This fact leads to the conclusion that 
the additional nonparabolicity of the hole dispersion law connected with the influence of the 
conduction band can be neglected for both narrow- and medium-gap semiconductor QDs. 

Taking into account the principal role of the dissymmetry coefficient C|, one can evaluate 
the influence of this parameter on the observable effects. With this purpose we calculate 
the lowest electron-hole pair energies as a function of for all QDHs under consideration 
(see Fig. |T0J). It is seen from Fig. |H] that when the parameter changes from —2 to 2, the 
corresponding energy differences E e _h(c^ = 2) — _E e _^(c^ = —2) constitute —175 meV for 
HgS/CdS (a = 2 nm), -15 meV for InAs/GaAs (a = 4 nm), and 20 meV for GaAs/AlAs (a = 

3 nm) QDHs. These differences should be quite accessible for the experimental detection. 
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VI. CONCLUSIONS 



The exact nonsymmetrized 8-band effective-mass Hamiltonian for an arbitrary 3- 
dimensional heterostructure has been obtained using the Burt's envelope-function represen- 
tation. The 2x2 electron and 6x6 hole energy-dependent Hamiltonians have been deduced. 
Within the spherical approximation, the 8x8, 2x2, and 6x6 radial Hamiltonians and the 
necessary BCs have been derived for spherical QDHs. The boundary conditions for radial 
symmetrized and nonsymmetrized Hamiltonians are different and lead therefore to different 
energy levels and wave functions. We have shown, further, that the CRs, which are com- 
monly used to match the solutions of the appropriate bulk kp Hamiltonians, coincide with 
BCs for the symmetrized Hamiltonians. A theoretical estimate for the value of the spin- 
orbit splitting of electron levels has been found. The energy levels of the nonsymmetrized 
8-band Hamiltonian have been calculated as a function of the dot radius for three spherical 
QDHs: a zero-gap semiconductor embedded into a wide-gap semiconductor (HgS/CdS), a 
narrow-gap semiconductor embedded into a medium-gap semiconductor (InAs/GaAs), and 
a medium-gap semiconductor embedded into a wide-gap semiconductor (GaAs/AlAs). It 
has been demonstrated that parameters of dissymmetry x( r ) an d £( r ), giving nonzero con- 
tribution to the multiband Hamiltonians only at the heterointerfaces, have, nevertheless, 
a strong effect on the electron and hole spectra. Thus, for practically important cases of 
relatively small QDHs with noticeably different effective-mass parameters of the constituent 
materials, the use of the obtained Hamiltonian is necessary for the adequate description of 
experiment. 
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APPENDIX A: INTERBAND MOMENTUM MATRIX ELEMENTS 

Using the general effective mass equations [Eq. (6.3) of Ref. [19| , one can see the origin 
of the interband momentum matrix elements t>i(r) and t>2(r) from Eq. ([!]): 

Aj Aj 

Vl (E,v) = -jr(S\p z \Z) --^(E-H^r'Hs^v) (u\p z \Z); (Al) 



M 

v 2 (E, r) = -- J2(E - H vv {v)Y l (S\p z \v) H vZ {v). (A2) 
n v 



To obtain Eqs. (|A1| ) and (|A2| ) it should be taken into account that within the developed 8- 



band approach the conduction band with the Bloch function \S) and the valence band with 
the Bloch functions \X), \Y), and \Z) are included explicitly, while all other bands with the 



Bloch functions \v) are considered to be remote. Further, following the technique of Ref. |19 



it is necessary to exclude the energy dependence from Eqs. (DO) and (|A2|) by replacing E 
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with an average energy, for instance the energy at the middle of the narrowest gap for the 
heterostructure compounds. Parameters t>i(r) and V2{r) are approximately considered to be 
constant in each layer of a heterostructure. 

In bulk, if Burt's material-independent basis functions coincide with the bulk Bloch 
functions, the second term in the r.h.s. of Eq. ( |A1| ) and the r.h.s. of Eq. ( |A2| ) vanish, because 
in this case Hs u (r) = and H v z{v) = 0. Therefore, one obtains v i = 2v {v = —2i(S\p z \Z)/h) 
and i>2 = 0, what results in £ = v (see Eq. (g)) and = 1 (see Eq. (g)). When materials 
constituting the heterostructure have close parameters, the second term in the r.h.s. of 
Eq. ( |A1[ ) and the r.h.s. of Eq. ( [^ ) are small compared with the first term in the r.h.s. of 
Eq. ( |A1| ). In this case £ does not differ significantly from v, and is close to 1. In a general 
case of disparate materials, can take arbitrary values. 

APPENDIX B: ENERGY-DEPENDENT SEPARATE HAMILTONIANS FOR 

ELECTRONS AND HOLES 

For narrow-gap semiconductors, the accurate way to take into account the coupling of 
conduction and valence bands, is to consider the 8-band Hamiltonian. However, sometimes 
it is easier to solve a CB or VB Schrodinger equation with energy-dependent effective- 
mass parameters. Solutions of these equations are just an approximation to the results of 
the 8-band model. In what follows, we deduce the 2x2 energy-dependent Hamiltonian for 
an electron and the 6x6 energy-dependent Hamiltonian for a hole from the exact 8-band 
nonsymmetrized effective mass Hamiltonian. 



1. 2x2 energy-dependent Hamiltonian for an electron 

a. Nonsymmetrized CB Hamiltonian 

We start with the nonsymmetrized 8-band Hamiltonian H defined by Eq. ([12]). The 

wave function i.e. a vector of eight envelope functions . . . , \l/8 is an eigenfunction of 
the matrix Schrodinger equation 

Hm = Em, (Bl) 

where E is an eigenenergy. To find the CB Hamiltonian, one should treat all VBs as remote. 
Therefore, we should exclude all VB envelopes, i.e. \l/3, . . . , \l/8, from Eq. (|B~1|). As seen, 
this exclusion is possible only within the approximation 71 = 0, 72 = 0, and 73 = 0, in other 
words, when the contributions to the hole effective-mass parameters from the remote bands 
(all bands except two CBs and six VBs) is negligible. This is a very close approximation, 
because the parameters 71, 72, and 73 are small for almost all materials and, determining 
contributions to the VB, they certainly have small influence on the electron levels. Under 
this approximation x — ~ 1/3 (see Eq. (|r|)) and it cancels from the Hamiltonian. Another 
necessary approximation is eg = —1, i.e. £ = —v, and therefore Vi = 0, v 2 = 2v (see Eqs. (^) 
and (@)). This is the only approximation that does not lead to the discontinuity of CB 
envelopes and ^2 at the heterointerface. Now we can express six VB envelopes ^3, . . . , 
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\l/8 in terms of two CB envelopes and ^2 using the last six equations of the set ( |B 1|) . 
This procedure results in 




3 



V 1 r 

^6 = — /c+^2, 

E 6 V 

* 8 = ^7*(V! Ml -^ A 



(B2) 



where e = 2m E/h 2 . Substituting the envelopes ( |B~2"D into the first two equations of the 
set (pl|), one obtains the sought CB Hamiltonian for the electron envelopes \&i and ^2- This 
Hamiltonian has the form 

Pe = ^ + (^) + L + L (5) - * c(e) ) ^+ + i? 

C e = -V2 (fc* g c {e) k- - L g c {e) k) , (B4) 

m v 2 ( 2 

a + — h 



where 



m c (e) 3 \e — e v e — e v + 5 



v 2 



Here, m /m c (e) is the inverse of the energy-dependent effective mass of an electron, and g c {e) 
is an energy- dependent interfacial parameter, which vanishes when the spin-orbit splitting 
5 is zero. Like x? this parameter gives a nonzero contribution to the Hamiltonian only 
at the heterointerface. The parameter g c {s) is responsible for the nonsymmetrical form of 
the Hamiltonian ( B3|) and for the mixing of the envelopes ^1 and ^2- When one solves the 
Schrodinger equation for an electron using the energy-dependent Hamiltonian ([B3|), one finds 
the eigenenergy E and eigenfunctions and ^2- Then, substituting these eigenfunctions 
into Eq. (B"2"|), one obtains the rest envelopes ^3, and therefore While the 

functions ^1 and ^2 are continuous at the heterointerface, the functions ^3, . . . , \l/8, as the 
envelope functions of all the other remote bands, are not.0 Note, that only the wave functions 
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i.e. vectors of eight envelope functions . . . , ar e orthonormalized. The envelopes 
and ^2 axe neither orthogonal nor properly normalized. When the nonparabolicity is 
not strong, in other words, when m c (e) only weakly depends on the energy, it is possible 
to choose one appropriate value of the energy, e. g. e , to find m c (e ) and g c (s ), and to 
substitute them into the Hamiltonian (|B3|). In such a way one obtains the Hamiltonian 
H e (eo). The eigenfunctions \&i and ^2 of this Hamiltonian will be orthonormalized and the 
rest six envelopes will be no longer needed. 



b. CB Hamiltonian and BCs for a spherical QDH 
The radial CB Hamiltonian nj?} for spherical QDHs can be derived from the Hamiltonian 



by the same way as the radial Hamiltonian ( ]2j ) has been obtained from the Hamiltonian 
(|12|) in the spherical approximation (see Sec. III). Thus we find 

«8 = (B6) 



where 



VS = ^ ^ ^ - A!" (M) (B7) 

and the operator A^\f3) is defined by Eq. (p6|) . Inside the z-th spherical layer, the Hamil- 
tonian ( PB| ) takes the form 

where A; is the spherical Laplacian and m c ,i(£) is the energy-dependent CB mass of the z-th 
material. Further, one should solve the Schrodinger equation with the Hamiltonian ( |B8| ) 
for each spherical layer and match the obtained solutions at the spherical heterointerfaces 
using the BCs ( ffip (see Sec. IV). The radial component of the CB current operator is 
obtained from the Hamiltonian ( P6| ) by the same way as the radial component of the current 
operator has been obtained from the Hamiltonian (E5). Thus, 




° - p(j + 1/2 -p) 9 -^]. (B9) 



In a two-layer spherical QDH the electron energy depends on the difference y C) i(e)— g c ,2{^) 
(as seen from the BCs (j3iP, (p9|)), where the indices "1" and "2" denote the interior and 
exterior materials, correspondingly. The value of this difference is usually very small for 
typical QDHs. Therefore, in the first approximation, one can find the energy spectrum E\ 
(I = j —p/2) neglecting the term proportional to p(j + 1/2 — p) (^g c ,i(s) — g c ,2(£)) in the BCs. 

Then, including this term as a perturbation, one finds the energy spectrum E{. It is seen 
that the energy levels with I = remain unchanged, while each energy level E\ with I > 1 
splits into two levels: E\ +l ^ 2 and E\ 1 ^ 2 . For the electron levels that are not very close to 
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the CB minimum (it is the case for the QDHs under consideration), the following estimate 
can be obtained 



(/ + 1) E l , +1/2 + 1 E l ~ 1/2 , +1/2 ,_ 1/2 h 2 (2l-l) f . . , .n . 



where = 2m Ei/fi . In Table [III] we have used Eq. ( |B10[ ) to estimate the spin-orbit 



splitting of the lowest P and .D levels. It is seen that the value of the splitting of the 
electron levels is of the order of 3 meV in all considered QDHs, and therefore this splitting 
can be neglected. This fact does not imply that the dependence of the parameters of the 
CB Hamiltonian on the energy can be neglected, too. As seen from Table [TIT], the electron 
effective masses in QDHs can differ by a factor of 2 from their values in the corresponding 
bulk materials. 



2. 6x6 energy-dependent Hamiltonian for a hole 

a. Nonsymmetrized VB Hamiltonian 

The deduction of the VB Hamiltonian is analogous to the deduction of the CB Hamilto- 
nian with the only difference: one should treat two CBs as remote. In order to express the 
CB envelopes and ^2 i n terms of the VB envelopes ^3, . . . , and to exclude them from 
the Schrodinger equation (|BT|), one should apply the approximation a = and eg = 1. This 
approximation has the same grounds as the approximation used above to obtain the CB 
Hamiltonian. Now, we express the CB envelopes ^1 and ^2 in terms of the VB envelopes 
\l/ 3 , . . . , from the first two equations of the set ( [BID and substitute them into the last 
six equations of the same set. As a result we have the VB Hamiltonian H^, which coincides 



with the Hamiltonian H (see Eq. fll2 )) where the first two rows and the first two columns 



are deleted and the effective-mass parameters are changed in the following way: 

7i li( £ ) = 7i + 



v 2 



3{e c - s) 



7,3 72,3 + ^Ty (BID 

Here, 7$(e) are the energy-dependent Luttinger parameters. In conformity with Eq. (0), 
one should change the parameter of dissymmetry x as follows 

X^X L {e) = X + f (B12) 
6{e c - e) 

As a result of the change ( |B12j ), the parameter of dissymmetry increases (see Table [TTT]), 
therefore the results of the symmetrized Hamiltonian (with \ L — 0) will deviate sharply from 
the exact solutions. The parameters 7 4 L usually weakly depend on the energy. Consequently, 
to obtain the hole spectrum one can use the Hamiltonian H^eo), where Eq is an average 
hole energy. 
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b. VB Hamiltonian and BCs for a spherical QDH 



The radial VB Hamiltonian wjfj for a spherical QDH coincides with the radial Hamilto- 
nian ( pip in which 7x — > 7f(e) and x - > (i n conformity with Eqs. ( Bll| ) and (pi2| )), 

v 2 

7 — > 7 L (e) = 7 H t r (in conformity with Eqs. ( |B11| ) and (0)) and where the first 

two rows and the first two columns are deleted. For the radial components of the hole wave 
function one should use the BCs (BTj), in which the radial component of the current operator 
Slfj is given by the matrix ( pl2"D where the first row and the first column are deleted and the 
parameters 71, 7, x are replaced by the parameters 7^ , 7 L , x L 1 correspondingly. 
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TABLES 



TABLE I. The 8-band effective-mass parameters of some III-V materials. The VB offset E v is 
chosen to be zero in GaAs. The parameters a, 71, 7, and x °f the spherical model are calculated 
from the listed here effective-mass parameters. 



Parameters 


GaAs 


AlAs 


InAs 


m c (m ) 


0.0665 c 


0.150 c 


0.02226 




7.10 d 


3.76 d 


19.67 c 




2.02 d 


0.90 d 


8.37 c 


7 3 L 


2.91 d 


1.42 d 


9.29 c 


E p (eV) 


28.0 b 


21. l c 


22. T 


E g (eV) 


1.519 c 


3.130 c 


0.418 c 


A(eV) 


0.341 c 


0.275 c 


0.380 c 


E v (eV) 





-0.532 d 


0.186° 


a 


-2.27 


0.11 


0.24 


7i 


0.96 


1.51 


1.97 


7 


-0.52 


0.09 


0.07 


X 


-1.52 


-0.69 


-0.87 



a Ref. g? 
b Ref. H 
c Ref. |9 
d Ref. 
e Ref. (13 
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TABLE II. The 8-band effective-mass parameters of some IV- VI materials. The VB offset E v 
is chosen to be zero in HgS. The parameters of the spherical model a, 71, 7, and x f° r CdS and 
X for HgS are calculated from the listed here effective-mass parameters. The parameters m c , 7^, 
and 7^ are not presented for HgS, because in a semimetal the band structure is inverted and these 
parameters do not have their original sense. 



Parameters 


HgS 


CdS 


m c (m ) 




0.18 f 


-y£ 

Q 




1.71 f 


1 L 




0.62 f 


E p (eV) 


13.2 a 


21. e 


E g (eV) 


-0.190 a 


2.56 c 


A(eV) 


0.07 c 


0.07 c 


E v (eV) 





-0.93 d 


a 


-1.0 C 


-2.57 


7i 


0.35 b 


-1.02 


7 


-0.67 b 


-0.75 


X 


-1.57 


-1.24 



a Ref. 
b Ref. H 
c Ref. H 
d Ref. [45 
e Ref. |6 
f Ref. 
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TABLE III. The spin-orbital splitting of electron energy levels (c^ = —1); the energy-dependent 
electron effective masses (eg = — 1); the energy shift of electron and hole levels due to finite values 
of x (at c i = 0)i an d the difference of \ L i n adjacent materials for the 6x6 model (at C£ = 1). 
If not indicated explicitly, / 0. E e and are the electron and hole ground state energies 
corresponding to the states lS 1 ^ ana - 1^3/2 • indices "1" and "2" denote the interior and 

exterior materials, correspondingly. 





HgS/CdS 


InAs/GaAs 


GaAs/AlAs 




a = 2 nm 


a = 4 nm 


a = 3 nm 


E (e) - E (e) (meV) a 

3/2 lr l/2 

E (e) - £ ( e) (meV) a 

5/2 ^3/2 


2.6 (2.0) 


3.0 (1.4) 


2.7 (3.0) 


2.0 (2.3) 


3.8 (2.2) 


4.5 (7.1) 


m c ,x(E e ) (m ) b 


0.056 ( - ) 


0.038 (0.022) 


0.083 (0.067) 


m C) 2(E e ) (m ) b 


0.097 (0.180) 


0.040 (0.067) 


0.115 (0.150) 


" ^e| x=0 (meV) c 


16.1 (6.3) 


-3.5 (-3.1) 


4.8 (7.0) 


^ " ^l x=0 ( m eV) c 


9.3 (6.3) 


-4.1 (-3.1) 


4.1 (7.0) 


X [{E h )-xk{Eh) d 


36.28 (-0.33) 


4.86 (0.65) 


0.73 (-0.83) 



a The theoretical estimate based on the 2x2 energy-dependent Hamiltonian for an electron (see 
Appendix B) is given in parentheses. 

b The corresponding bulk effective mass (see Tables | and O) is given in parentheses. 

n 2 

The result of the empirical estimate E e r h \ — E e i h \ = n(Xl ~ X2) is given in parentheses. 

y ' y ' x=o moa 2 

d The difference xi ~ X2 for the 8-band model (see Tables | and |0|) is given in parentheses. 



c 
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FIGURES 




(fjgS)'- 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 ' 1 1 1 1 



1.0 1.5 2.0 2.5 3.0 

a (nm) 

FIG. 1. All the discrete electron energy levels in the HgS/CdS QDH as a function of the 
quantum dot radius. The and D^L energy levels are not shown here and in Figs. ^ and || 

because in the chosen scale they coincide with the levels P^X an d ^3/2 > correspondingly. Solid lines 
represent the result of the symmetrized 8-band model (ce = 0, x = 0). With the nonsymmetrized 
valence band part of the Hamiltonian (x 7^ 0), dashed lines show the case eg = while gray bands 
represent the continuous change of from 1 to —1. The gray bands refer to a possible variation in 
energy due to conduction band/valence band coupling via the position dependence of the interband 
momentum matrix element. Here and in Figs. |||-|^, the insert with inscriptions = — 1 and C£ = 1 
shows to which values of the edges of the gray bands are related. 
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FIG. 2. All the discrete electron energy levels in the InAs/GaAs QDH as a function of the 
quantum dot radius. Other denotations are the same as in Fig. [j]. 
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FIG. 3. All the discrete electron energy levels in the GaAs/AlAs QDH as a function of the 
quantum dot radius. Other denotations are the same as in Fig. 
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FIG. 5. All the discrete hole energy levels with j = 3/2 in the InAs/GaAs QDH as a function 
of the quantum dot radius. Other denotations are the same as in Fig. [j]. 
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a (nm) 

FIG. 6. All the discrete hole energy levels with j = 3/2 in the GaAs/AlAs QDH as a function 
of the quantum dot radius. Other denotations are the same as in Fig. |l]. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

r (nm) 

FIG. 7. S-type radial components of the wave functions of the electron and hole ground states 
in the HgS/CdS QDH (radius a = 2 nm) within the nonsymmetrized 8-band model. Solid and 
dashed lines denote the cases = 1 and = —1, correspondingly, while gray bands represent the 
continuous change of eg within these limits. Each radial wave function is normalized by unity, i.e. 
the integral probability J2fi Io° r 2 R^(r)dr = 1, where \i labels the radial components. Contributions 
to the integral probability from the depicted radial components vary from 68.6 % to 64.8 % for 
an electron and from 74.3 % to 81.1 % for a hole when eg changes from 1 to —1. At the same 
time the electron energy changes from 533.6 meV to 678.8 meV and the hole energy changes from 
-247.2 meV to -202.1 meV. 
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1 2 3 4 5 6 

r (nm) 

FIG. 8. S-type radial components of the wave functions of the electron and hole ground states in 
the InAs/GaAs QDH (radius a = 4 nm) within the nonsymmetrized 8-band model. Normalization 
of each radial wave function and denotations are same as in Fig. Contributions to the integral 
probability from the depicted radial components vary from 80.6 % to 79.2 % for an electron and 
from 75.3 % to 78.7 % for a hole when changes from 1 to —1. At the same time the electron 
energy changes from 736.0 meV to 757.1 meV and the hole energy changes from —73.6 meV to 
-68.2 meV. 



32 




12 3 4 

r (nm) 

FIG. 9. S-type radial components of the wave functions of the electron and hole ground states in 
the GaAs/AlAs QDH (radius a = 3 nm) within the nonsymmetrized 8-band model. Normalization 
of each radial wave function and denotations are the same as in Fig. [7[ Contributions to the integral 
probability from the depicted radial components vary from 85.4 % to 86.9 % for an electron and 
from 88.9 % to 86.6 % for a hole when changes from 1 to —1. At the same time the electron 
energy changes from 1880.4 meV to 1850.4 meV and the hole energy changes from —111.6 meV to 
-122.9 meV. 
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The lowest electron- hole pair energies in different QDHs as a function of eg (x / 0). 
related to the closest curve and indicates the corresponding result of the symmetrized 

: 0, X = 0). 



34 



REFERENCES 



t Electronic address: pokatilov@add.moldova.su 

* Electronic address: fonobero@usm.md 

§ Permanent address: Laboratory of Multilayer Structure Physics, Department of Theoret- 
ical Physics, State University of Moldova, A. Mateevici 60, MD-2009 Chi§inau, Moldova. 
Electronic address: fomin@uia.ua.ac.be 

* Also at: Universiteit Antwerpen (RUCA), Groenenborgerlaan 171, B-2020 Antwerpen, 
Belgie and Technische Universiteit Eindhoven, P. O. Box 513, 5600 MB Eindhoven, The 
Netherlands. Electronic address: devreese@uia.ua.ac.be 

1 J. M. Luttinger and W. Kohn, Phys. Rev. 97, 869 (1955). 

2 A. Baldereschi and N. O. Lipari, Phys. Rev. B 8, 2697 (1973). 

3 J. B. Xia, Phys. Rev. B 40, 8500 (1989). 

4 V. M. Fomin, V. N. Gladilin, J. T. Devreese, E. P. Pokatilov, S. N. Balaban, and S. N. 
Klimin, Phys. Rev. B 57, 2415 (1998). 

5 G. B. Grigoryan, E. Kazaryan, Al. L. Efros, and T. V. Yazeva, Sov. Phys. Solid State 32, 
1031 (1990). 

6 A. I. Ekimov, F. Hache, M. C. Schanne-Klein, D. Ricard, C. Flitzanis, I. A. Kudryavtzev, 
T. V. Yazeva, A. V. Rodina, and Al. L. Efros, J. Opt. Soc. Am. B 10, 100 (1993). 

7 C. R. Pidgeon and R. N. Brown, Phys. Rev. 146, 575 (1966). 

8 P. C. Sercel and K. J. Vahala, Phys. Rev. B 42, 3690 (1990). 

9 Al. L. Efros and M. Rosen, Phys. Rev. B 58, 7120 (1998). 

10 F. M. Gashimzade, A. M. Babaev, and M. A. Bagirov, J. Phys.: Condens. Matter 12, 
7923 (2000). 

11 P. C. Sercel, Al. L. Efros, and M. Rosen, Phys. Rev. Lett. 83, 2394 (1999). 

12 C. Pryor, Phys. Rev. B 57, 7190 (1998). 

13 O. Stier, M. Grundmann, and D. Bimberg, Phys. Rev. B 59, 5688 (1999). 

14 J. M. Shi, V. N. Freire, and G. A. Farias, Eur. Phys. J. B 14, 337 (2000). 

15 Ch.-Y. Hsieh, D.-S. Chuu, and Ch.-Y. Hsieh, J. Phys.: Condens. Matter 12, 8641 (2000). 

16 W. Jaskolski and G. W. Bryant, Phys. Rev. B 57, R4237 (1998). 

17 E. P. Pokatilov, V. A. Fonoberov, V. M. Fomin, S. N. Klimin, and J. T. Devreese, Bullet. 
Amer. Phys. Soc. 45, 988 (2000). 

18 M. V. Kisin, B. L. Gelmont, and S. Luryi, Phys. Rev. B 58, 4605 (1998). 

19 M. G. Burt, J. Phys.: Condens. Matter 4, 6651 (1992). 
20 M. G. Burt, Phys. Rev. B 50, 7518 (1994). 

21 M. G. Burt, J. Phys.: Condens. Matter 11, R53 (1999). 

22 B. A. Foreman, Phys. Rev. Lett. 80, 3823 (1998). 

23 G. Bastard, Phys. Rev. B 24, 5693 (1981). 

24 Y. R. Lin-Liu and L. J. Sham, Phys. Rev. B 32, 5561 (1985). 

25 G. A. Baraff and D. Gershoni, Phys. Rev. B 43, 4011 (1991). 

26 C. Y.-P. Chao and S. L. Chuang, Phys. Rev. B 46, 4110 (1992). 

27 G. Edwards, E. C. Valadares, and F. W. Sheard, Phys. Rev. B 50, 8493 (1994). 

28 B. A. Foreman, Phys. Rev. B 48, 4964 (1993). 

29 B. A. Foreman, Phys. Rev. B 56, R12748 (1997). 

30 R. van Dalen and P. N. Stavrinou, Semicond. Sci. Technol. 13, 11 (1998). 
31 F. Mireles and S. E. Ulloa, Phys. Rev. B 60, 13659 (1999). 



35 



32 F. Mireles and S. E. Ulloa, Phys. Rev. B 62, 2562 (2000). 

33 A. T. Meney, B. Gonul, and E. P. O'Reilly, Phys. Rev. B 50, 10893 (1994). 

34 S. De Franceschi, J. M. Jancu, and F. Beltram, Phys. Rev. B 59, 9691 (1999). 

35 L. W. Wang, Phys. Rev. B 61, 7241 (2000). 

36 D. J. BenDaniel and C. B. Duke, Phys. Rev. 152, 683 (1966). 

37 S. Adachi, Physical Properties of III- V Semiconductor Compounds (John Wiley & Sons, 
New York, 1992). 

38 S. L. Chuang, Physics of Optoelectronic Devices (John Wiley & Sons, New York, 1995). 

39 Landolt-Bornstein, Numerical Data and Functional Relationships in Science and Technol- 
ogy, New Series (Springer, Berlin, 1982), Vol. III/17a. 

40 p r0 p er ti es of Alluminium Gallium Arsenide, edited by S. Adachi (INSPEC, London, 1993). 

41 K. Dybko, W. Szuszkiewicz, E. Dynowska, W. Paszkowicz, and B. Witkowska, Physica B 
256, 629 (1998). 

42 71 + 27 = —1.0 is taken from Ref. O and 71 — 27 = 1.696 is taken from Ref. H3. 



43 M. von Truchsess, A. Pfeuffer-Jeschke, C. R. Becker, G. Landwehr, and E. Batke, Phys. 
Rev. B 61, 1666 (2000). 



44 P. Lawaetz, Phys. Rev. B 4, 3460 (1971). 

45 A. H. Nethercot, Phys. Rev. Lett. 33, 1088 (1974). 



46 C. Hermann and C. Weisbuch, Phys. Rev. B 15, 823 (1977). 

47 Parameters m c , 7^, and 7 L are computed applying the quasicubic model@ (i.e. l/m c = 
(2/m c± + l/m c ||)/3, ^ = (2/mj_ + l/mf)/3, 7 L = {l/mj_ - l/mf)/3) to the effective 
masses of wurtzite CdS. m c j_, m c \\, mj_, mj 4 were determined in Ref. |49|, using the empirical 
pseudopotential method. 

48 Yu. M. Sirenko, J. B. Jeon, K. W. Kim, M. A. Littlejohn, and M. A. Stroscio, Phys. Rev. 
B 53, 1997 (1996). 

49 J. B. Xia and J. Li, Phys. Rev. B 60, 11540 (1999). 



36 



